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Abstract 

We investigate non-Gaussianity in general single field models without assuming 
slow-roll conditions or the exact scale-invariancc of the scalar power spectrum. The 
models considered include general single field inflation (e.g. DBI and canonical in- 
flation) as well as bimetric models. We compute the full non-Gaussian amplitude A, 
its size /nl; its shape, and the running with scale ung- In doing so we show that 
observational constraints allow significant violations of slow roll conditions and we 
derive explicit bounds on slow-roll parameters for fast-roll single field scenarios. A 
variety of new observational signatures is found for models respecting these bounds. 
We also explicitly construct concrete model implementations giving rise to this new 
phenomenology. 
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1 Introduction 

Non-Gaussianity as probed by upcoming experiments such as Planck will be a very powerful 
tool for selecting models of primordial structure formation [U [21 [21 [H [3, IS]- One of the 
most interesting questions one may hope to answer is: Are single fleld models of structure 
formation sufficient, or do we need multiple fields to meet observational constraints? The 
hallmark of single field models that generate large non-Gaussianity is often thought to be an 
equilateral type non-Gaussian shape [HIHl IHl [IDl [HI 112] ■ Known modifications to this shape 
can be achieved by introducing, for example, non-Bunch-Davies vacua [7], [13], [TU [TJl US] 
or resonant running (ISl [E] [iHl [IHl [20] • However, models with additional fields, such as 
muhifield infiation [2Il[221[23l[2a[25l[26l[2Zl[2Hl[2ai3Ql[3ll|321[Sll3ll|3^ models adding 
curvaton fields [361 [33 I3H1 [311 [IS], HH [121 [131 [11] , can give rise to a much greater variety 
of non-Gaussian signatures (see e.g. [33])- It is therefore interesting to investigate whether 
multifield models are really unique in predicting all such alternative shapes. 
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In this paper we demonstrate that non-slow-roll single field models can produce a much 
richer variety of detectable non-Gaussian signatures than those following from slow-roll condi- 
tions. Typically slow-roll conditions are expected to be imposed by observational constraints 
derived from primordial structure formation [1] and from bounds arising by requiring solu- 
tions of e.g. fiatness, homogeneity and isotropy problems. In [8j it was shown that for exactly 
scale-invariant non-canonical models, a much wider range of slow-roll parameters is allowed 
than those satisfying slow-roll conditions, whilst still satisfying all observational bounds. We 
extend this work by dropping the requirement of exact scale- invar iance, explicitly comput- 
ing the corresponding new bounds for slow-roll parameters and showing that an even wider 
range of parameters becomes observationally viable. 

Non-Gaussianities for non-slow-roll single field models have been previously computed 
for the case of exact scale invariance = 1 |8] and that of an effective DBI type action [llj. 
We will generalize this work, restricting ourselves to expanding solutions as provided, for 
example, by inflationary and bimetric models [IHl E]. Contracting models have been 
considered in the context of ekpyrotic/cyclic models [HI HE] and effective perfect fluid mod- 
els [ISl |50] .[] In this paper we hope to complete the picture in the expanding case by deriving 
non-Gaussianities for general single fleld models without assuming smallness of slow-roll pa- 
rameters ej or exact scale invariance (n^ = 1). We investigate the phenomenology associated 
with such models, showing that they have several distinct observable signatures. These 
include a generically large running of non-Gaussianity ung with scale, the suppression of 
the size of non-Gaussianity /nl by non-slow-roll and the emergence of qualitatively new 
shapes for single fleld models e.g. ones with dominant contributions in the folded limit 
2ki = 2k2 = k^. These considerations widen the number of potential future observations 
that could be explained by single fleld models. And even if future measurements turn out 
to favor multifleld models, it is still important to fully understand the potential effects of 
individual flelds that may participate in such models. 

The plan of the paper is as follows. In section [2] we set the stage by considering the two- 
point correlation functions for scalar and tensor perturbations in general single fleld models. 
We compute the resulting observational constraints on slow-roll parameters and also summa- 
rize bounds coming from other sources. Then, in section |3| we compute non-Gaussianities, 
explicitly giving the full amplitudes for general single fleld models without assuming slow-roll 



conditions or exact scale-invariance (section [3.1 ). Observational signatures are considered 



in detail for /nl (section 3.2[ ), non-Gaussianity shapes (section 3.3) and running with scale 



f^NG (section 3.4). In section |4] concrete models are provided for the phenomenologically in- 
teresting cases by explicitly constructing their actions. In section [5] we conclude, summarize 
results and point out intriguing possibilities for further work. 



^While this paper was being finished, references [51] [52] appeared, which further investigate contracting 
models. |52| also considers strong coupling constraints for expanding models. 
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2 Single field models without slow-roll 



The models we shall consider have an action of the form [53t IM] 

'R 



S 



d x\/—g 



(2.1) 



where the pressure P is a general function of a single scalar field and its kinetic term 
X = —\g^"d^j_(j)dy(l). For a canonical field P[X,(j)) = X — Vi^cj)). In principle P can contain 
higher derivatives of the field as well, however such terms are expected to be suppressed 
by the UV cut-off scale in an effective field theory. This motivates studying an action of 



the form (2.1), which is the most general Lorentz invariant action for a single scalar field 



minimally coupled to gravity that contains at most first derivatives of the field. The pressure 
p and energy density p for these models are given by 



P = Pix,<p), 

p = 2XP^x~P{X, 



(2.2) 



where = dP/dX. 

We will find it useful to characterize the non-canonical nature of (2.1 ) with the following 
quantities |55] 

2 P.x 



P X + 2XPxx 



XP,x + 2X^P^xx 



X^P^xx + -^X^Pxxx , 



(2.3) 
(2.4) 
(2.5) 



where Cg is the adiabatic speed of sound. An infinite hierarchy of slow-roll parameters e and 



and their derivatives can then be constructed for a theory (2.1) with an FRW metric with 
scale factor a{t) and corresponding Hubble rate H{t) = a/a. These hierarchies are 



H 

IP' 



V 



e 

7h' 



Vs 



(2.6) 



The slow-roll approximation requires that e,rj...es,ris... <^ 1 and accelerated expansion takes 
place when e < 1. Note that these definitions are more general than alternative definitions 
in terms of the scalar potential and its derivatives, which assume canonical kinetic terms. 



Perhaps more fittingly (2.6) have therefore also been dubbed "slow- variation parameters" 



ETJEH]- DBI inflation |59], for example, is slow varying in this sense yet allows steep potentials 
and hence the infiaton is not slow-rolling. However, in order to agree with convention and 



avoid confusion, we will refer to (2.6) as slow-roll parameters from now on. 



Truncating slow-roll parameters at some derivative order, setting the remainder to zero, 
allows us to approximate H and Cg with Taylor expansions. Following [HI E] we consider 



models where slow-roll conditions are broken by e,es ~ We will assume that e,es are 

approximately constant, so that slow-roll conditions still hold for the remaining negligible 
parameters rj^rjs... ~ 0. 



2.1 Observational constraints on (non)-slow-roll 



Large classes of both inflationary and bimetric models can be described by action (2.1). 



Bimetric models generically have e > 1 and are discussed in detail in section |4.3 For the 



inflationary branch e < 1, however, observations can be used to constrain the allowed values 
of slow-roll parameters e, e^. These constraints are derived throughout the paper and we 
summarize the resulting bounds upfront here. The observational bounds considered are the 
following: 

• The spectral index n^: The spectral index of scalar perturbations is computed 
in section [2I2 In terms of the slow-roll parameters it can be written as 



1 2e + e. 

The WMAP bounds on in the presence of tensor perturbations are Ug = 0.973±0.028 
at 2cr. Note that these constraints are dependent on e.g. the precise value of Hq 
and the physics of reionization [20] . Contrary to the canonical field case, for a general 



single field model (2.1) the smallness of ri^ — 1 ^ 1 no longer automatically translates 
into a constraint on e alone. Instead, for any parameter e, one can find a corresponding 
speed of sound profile parameterized by e^, which returns any desired spectral index. 
Bounds on Ug therefore only translate into constraints on the relation between e and 
e^, not on either variable by itself. 

The tensor-to-scalar ratio r: The tensor-to-scalar ratio r is similarly constrained 
by CMB observations. The WMAP experiment provides a bound r < 0.24 at 2a [1]. 



In section 2.3 we discuss how the first order relation r = 16eCs is modified in general 
single field models without slow-roll and compute corresponding bounds on slow-roll 
parameter e. We find 

e < 0.4. (2.8) 

The running ung- Introducing non-canonical kinetic terms can lead to a strong 
running of non-Gaussianities with scale rij^G ^ EB ESI IS3], irrespective of whether the 
2-point correlation function as measured by Ug is near-scale- invariant or not. This can 
impose further constraints on e. If non-Gaussianities show a blue tilt, this means they 
peak on smaller scales. For a perturbative treatment to be applicable, one therefore 
needs to ensure that non-Gaussianities on the smallest observable scales are not too 



large for a perturbative approach to break down. In section 3.4 we identify under 



which conditions non-Gaussianities are blue tilted and what constraints on e follow. 
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We show that constraints are in principle model-dependent. However, for a very large 
subclass of models one can obtain rough bounds 



e < 0.3 if /nl(CMB) ~ O(IOO) (2.9) 
e < 0.5 if /nl(CMB) ~ (2.10) 

This result agrees with bounds obtained in the exactly scale invariant case However, 
we show that these limits are in fact relatively insensitive to the precise value taken by 
Hs, but strongly depend on parameters c^, S, A. 

The Big Bang problems: Observational bounds only place constraints on primordial 
perturbations for a small range of scales: CMB (~ lO'^Mpc) to galactic scales (~ IMpc). 
ns,r and unc as discussed above are therefore only observationally constrained over 
these scales. 

In an inflationary setup this corresponds to approximately 10 e-folds, where the number 
of e-folds of inflation is given by dN = —Hdt. In order to resolve the flatness, 
homogeneity and isotropy problems of cosmology, a greater number of e-folds is needed. 
Conventionally this is chosen ~ 60, although the minimal amount of e-folds needed 
is not very well defined and depends e.g. on the energy scale of inflation, details of the 
reheating mechanism, etc. |6l] The inflationary condition e < 1 is sufficient to yield 
an attractor solution towards flatness, homogeneity and isotropy [8J. A large number 
of e-folds requires e < 1 for several Hubble times, which is controlled by parameter rj. 
As we consider cases for which ~ ^ 1, the fractional change in e per Hubble time 
is small and hence starting with e < 1 guarantees that a prolonged phase of inflation 
takes place. 

However, we should note that several other solutions are possible as well. For example, 
the scaling phase responsible for near-scale-invariant perturbations on observable scales 
may be followed by a different inflationary phase (scaling or not) [H]. Particle physics 
may suggest that inflation is actually a multiple step process, where one relatively 
short initial phase is responsible for observable structure on very large scales [65j . 
Inflation then ends and is followed by (an)other phase(s) of accelerated expansion 
later. It could also be that an altogether different mechanism[^ resolves the problems in 
question here [HSl EZl EH]- In any case no competitive constraints on e can be derived 
from consideration of these cosmological problems. Also note that bimetric models 
resolve these problems in a radically different way [66], ETj |68] . 



Inflationary models of the type envisaged by (2.1 ), with e < 0.3, meet all the observational 
bounds. 



^And of course, as the flatness, homogeneity and isotropy problems are essentially initial value problems, 
the correct initial values may be set by pre-inflationary physics. Whilst such a solution may be very non- 
definitive and as such less appealing, one should keep in mind that there is certainly no guarantee that all 
cosmological puzzles can be answered at energy scales we have (in-)direct access to. 
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2.2 The spectral index 



In this section we sketch the calculation of the primordial power spectrum of scalar pertur- 
bations. Scalar perturbations in the metric 



ds^ = (1 + 2<i>)dt^ - (1 - 2<i>)a^{t)-f,jdx'dx\ 
can be analyzed in terms of the the gauge invariant curvature perturbation ( 

5p + 3p 2p $ 

3{p + p) 3{p + p)H' 



(2.11) 



(2.12) 



where p and p are the background pressure and energy density respectively. Perturbing (2.1 ) 
to quadratic order in ( one obtains 



^2 



^ / d'xdr 



dC 
dr 



(2.13) 



where r is conformal time, and z is defined as 2; = ay2e/cs. 

As a consequence of considering non-canonical kinetic terms, we have introduced a dy- 
namically varying speed of sound. This means it is convenient to work in terms of the 
"sound-horizon" time dy = c^dr instead of r [8]. For the case relevant here, i.e. rj = rjs = 0, 
one finds 

+ Dag - 

In terms of y-time the quadratic action takes on the familiar form 



^2 



Mil 



d xdy q 



C - (VC)^ 



c.z 



a^f2e 



(2.15) 



where ' = d/dy. It is useful to list the behavior in proper time r and in y-time of some 
relevant quantities for later reference: 



a ~ (— y) '^+' 



— 6 



(2.16) 



Upon quantization the perturbations are expressed through creation and annihilation 
operators as follows. 



C(Z/, k) = Mfc(y)a(k) + ul{y)a\-\^) 



[a(k), at(k')] = (27r)=^53(k - k'). (2.17) 



In terms of the canonically-normalized scalar variable v = M-piqQ the equations of motion 
for the Fourier modes are 

/ n" \ 

Vk = Q. (2.18) 
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If e,r],es,r]s are constant, as is the case here, one can solve for q"/q exactly 

'"-^C^-iV (2.19) 



q y"^ \ 4 

The precise solution for Vk{y) depends on which initial vacuum is chosen. Typically this is 
the Bunch-Davis vacuum, which is defined to be the state annihilated by a(k) as t — t- — oo. 
With this vacuum choice the solution for Vk{y) is 



Vk{y) = ^V^Hl'\-ky). (2.20) 

where hI'^^ are Hankel functions of the first kind. However, we note that different adiabatic 
vacuum choices are possible in principle and also have distinctive observational signatures 
(particularly non-Gaussianities) [HI [H [ISl E] • For the modes defined in (2.17) one finds 

My) = f}^'^ vk{y) = -S^.\F-^H^"\-ky)- (2.21) 



Therefore the 2-point correlation function for ( is: 



(C(ki)C(k2)) = i2n)'6\k, + k^)^ , (2.22) 
where the expression for the power spectrum is 

Note that we have approximated the propagator Uk{y) assuming — 1 ^ 1 in accordance 
with observations here (for details see appendix A.l - a more general expression for can be 
found in [52] )• The bar symbol means that the corresponding quantity has to be evaluated, 
for each mode fc, at sound horizon exit, i.e. when y = k~^. The spectral index is finally 
given by 

n,-l = — — ^ = 3 - 2z/ = . 2.24 

d\nk e, + e - 1 ^ ^ 

Therefore exact scale invariance obtains when = 1, i.e. eg = — 2e. Constraints on Ug 

consequently translate into relational constraints on e and e^. Importantly this means no 

direct bound on either e or can be obtained in this way, because for any parameter choice 

e one can find an appropriate to yield a (near) scale- invariant solution. This is in contrast 

to the case of canonical fields, where the smallness of — 1 requires e to be small. 



2.3 The tensor-to-scalar ratio r 

In inflation [45j primordial quantum tensor fluctuations are sourced on sub-horizon scales 
before leaving the horizon and freezing out in analogy to scalar fluctuations. Present ob- 
servational upper bounds on the tensor-to-scalar ratio r therefore constrain the amount of 
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gravity waves that could have been generated in the early universe. Here we use bounds on r 
to constrain slow-roll parameter e. Note that bimetric models of the primordial epoch do not 
resolve the horizon problem for gravitational waves [16]. As no quantum tensor fluctuations 
are amplified to become classical outside the horizon one does not expect any significant 
primordial gravitational wave contribution for these models. Such models therefore predict 
r ^ 0. 



The action governing tensor modes is [69 



S 



a 
a 



drj 



(2.25) 



where the sum runs over the two distinct polarization modes and we have introduced v 
to distinguish between variables v and v associated with tensor and scalar perturbations 
respectively. The action is written in terms of a Mukhanov-Sasaki variable ?)k,n = ^^^^^^^^k,n; 
where a is the scale factor of the FRW metric. Upon quantization v can be expressed in 
terms of creation and annihilation operators as 



k,n' 



The equation of motion for tensor modes can then be written as 



Vk = 0. 



(2.26) 



(2.27) 



Accordingly we can now write down the power spectrum of tensor perturbations, where the 
sum again runs over polarization modes 



(2.28) 



which can be found to be [571 ES] 



Ph = 2^^ 



-3 



r(/i) 



where 



r(3/2) 



(1-e) 



1 - e 



V2H 



ttM 



pi 



k=aH 



1 

2' 



(2.29) 



(2.30) 



Importantly the tensor power spectrum is evaluated at horizon crossing for the tensor modes. 



For general theories with a non-canonical kinetic term (2.1) this horizon is different from the 



corresponding horizon for scalar perturbations. The spectral index of tensor perturbations 
Tit then is 

-2e 

nt = . (2.31) 



1 - e 



The tensor-to-scalar ratio r is given by 

'1 - e)2M-i 



3 



[es + e 



r(/i) 



r(i) 



2 - 2 



16c.(A;c)e(^) (2.32) 



where has to be evaluated at scalar horizon crossing and and if^ correspond to the 
Hubble factor at freeze-out for tensor and scalar modes respectively Note that the 

dependence on the ratio of H at different horizon crossings can induce a strong scale depen- 
dence for r. 

Present-day observational constraints limit r < 0.3. The WMAP 95 % confidence level 
bound is r < 0.24 For a given value of Ug this can therefore be used to constrain 

eCs{k(^). Ref |70] perform an MCMC sampling for general single field models ( 2.1[ ) to find 



eCs < 0.023 at 2a for a flat prior (note that the result is somewhat prior dependent). For 
a horizon crossing sound speed of Cg = 0.1 this corresponds to e < 0.23 which agrees with 
the approximate bound derived by |Bj. However, as we will show in the following section 
single field models with lower sound speeds Cg ~ 0.05 can satisfy current bounds on the size 
of non-Gaussianity, so that this bound on e can be relaxed to obtain 

e < 0.4 (2.33) 



for general single field models as considered here (2.1). 



3 Non-Gaussianity 

In this section we compute the 3-point correlation function for the curvature perturbation 
(. Perturbing (2.1) to cubic order in ( one obtains the effective action derived in [551 [7]. 



The result is valid outside of the slow-roll approximation and for any time-dependent sound 
speed: 



PI 



ae , 



S 1- 



1 



2A 



CI 



+ 



2e, + l 



, . c^)C(9C)^-2a-C(aC)(ax) 



3 + 3c^)CC' 



^ UiW « + 2-J'00x)S\ + ^JS'OOx? + 2/(0 



(3,i; 



As a comparison with ( 2.23 1 shows, we have assumed n 
^This bound is obtained assuming nt = 



1 = 3 — 2j^ ^ 1 here, i.e. near-scale-invariance. 
however. The relation between nt and r is of course more 



involved for general single field models as a comparison between (2.311 and (2.32) shows. However, given 



that we obtain a red tilt of tensor modes (and Ph therefore peaks on the largest scales) the bound on r 
should only weakly depend on the precise value of the tilt. 
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where dots denote derivatives with respect to proper time t, d is a spatial derivative, and x 
is defined as 



(3.2) 



Note that the final term in the cubic action, 2f{Qj^\i is proportional to the linearized 



5C 



equations of motion and can be absorbed by a field redefinition C Cn + /(Cn) - for details 



see appendix A.l 



At first order in perturbation theory and in the interaction picture, the 3-point function 
is given by [7T1 [55], |7j 



(C(t,kOC(t,k2)C(t,k3)) 



dt'([C(t,kOC(t,k2)C(t,k3),//int(t')]) 



(3.3) 



to 



where i^int is the Hamiltonian evaluated at third order in the perturbations and follows 
from (3.1). Vacuum expectation values are evaluated with respect to the interacting vacuum 

The 3-point correlation function is conventionally expressed through the amplitude A 



(C(ki)C(k2)C(k3)) = (27r)^(5^(ki + k2 + k,)P, 



1 



A. 



(3.4) 



Also by convention the power spectrum in the above formula is calculated for the mode 
K = ki+k2+k'i. By using (2.17) we can now calculate the amplitudes for each term appearing 
in the action (3.1). Using a self-evident notation the overall amplitude is therefore given by 

(3.5) 



A = Ari 



A, 



where accounts for the d^dxd^x {d'^Q^dxY terms in the action (3.1). 

In the following subsections we first present the main result (Sec. 3.1), which is a general 
expression for A where e is not constrained to satisfy e <^ 1 are given. We then study three 
physically interesting properties of these three-point functions: The size of the amplitude 
as parameterized by the variable fNL (Sec. 3.2), its shape (Sec. 3.3) and the running of the 
non-Gaussianity hng (Sec. 3.4). Current observational bounds are also discussed for these 
properties. 



3.1 The full amplitudes 

We first present results for the full non-Gaussian amplitude A. Details of the calculation 
are given in appendix A.l We stress that our results do not assume slow-roll conditions for 
e, and are therefore valid even when these parameters become e, ~ Providing the 

full amplitudes here will hopefully prove useful for the computation of non-slow-roll non- 
Gaussianities for specific single field models in the future. For the exactly scale-invariant 
case analogous non-slow-roll results were presented in [S] and for an effective DBI type action 
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these have been derived in fH] without assuming scale-invariance. However, as we will see 
in the following subsections, it is precisely the part of the amplitude that vanishes in DBI 
type models that is instrumental for much of the new phenomenology we find for general 
slow-roll violating single field models (2.1). 

The following variables will be useful in computing amplitudes 

2e + 2e — e.. 



ai = — 1 = 3 — 2z/ 

es + e-1 

Also note that we can write the parameter A as [551 E] 

'2fx - 1 



a2 



es + e-l 



(3.6) 



A 



"6 



where 



X dH 



X 



3e 



ex 



X 



(3.7) 



(3.^ 



and we note that there is in principle no dynamical requirement fixing ex (and hence fx) to 
be small or large, even in the presence of constraints for e [S^. For DBI models (discussed 
in 4.1) /^^^ = 1 "~" so consequently the first term in the action (3.1) vanishes and there 
is no contribution to the amplitude. For the computation of .4^3 we otherwise assume 
for simplicity that fx (and hence the kinetic part of e, ex) is constant. Computing the 
amplitude in this way for each term in the action (3.1) we obtain 

l)X^3(a2) + cl{e + es- l)I^s{ai 



A 



A 



C{9C)= 



A. 



A, 









[ 2K^ 




f kik2h 




[ 2K^ 




f kik2h 




[ 2K^ 




f kik2k2, 






1 


f kik2k 


16c2 


V 2K^ 



)ns-l 
[-e%c9x("2)] ; 



[{e + es-l)Ux 
[{e - 3)1^^2(02) + 3c^X^^2(ai 
[(e - 2e, + 1)1^(^9^2(02) - clX^(QQ2{ai 



(3.9) 



where 



2^/3 (a) 



an 

cos — r(3 + a 
2 



COS —1 (1 + a) 



(2 + «)^E^^'^J-(i + ")]^E^^'^: 
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%9C)2(«) 



cos 



i<j 



a — 



an 

■ cos -;5-r(l + a) 



3 



, ^ . A;iA;- - akxkikj, 



: cos -^r(l + a) 



avr. 



X,2(a) =cos— r(l + a)(2 + a/2) 



^k]~^kik]\1k^k2k^ 



(3.10) 



where we remind ourselves that A^^ accounts for the d(dxd'^x ^-^id {d'^C){dxY terms in the 



action (3.1 ). K = ki + k2 + k^ is the sum of the lengths of the three wavevectors ki, k2, ks. 



Note that, in deriving the above results, the propagator Uk{y) has been approximated in the 



\ky\ -C 1 limit (for details see appendix A.l, especially equation (A. 3)). The expressions 



derived here reduce to the correct amplitudes in the DBI [TT], scale-invariant |8] and slow- 
roll [7] limits. 



3.2 / 



NL 




0.05 0.10 0.15 0.20 0.25 



Figure 1: /nl is plotted for parameters Left: = 1, fx = 0.01 and Ug = 0.96, 1, 1.04 from 
top to bottom respectively. Middle: Cs = 0.05, fx = 0.01 and Ug = 0.96, 1, 1.04 from top to 
bottom respectively. Right: = 1, /x = 3 ■ 10^ and = 0.96, 1, 1.04 from top to bottom 
respectively. 

The size of the three-point correlation function is usually quoted using parameter /nl- 
In terms of the amplitudes A this is conventionally defined as [8] : 



(3.11) 



where we have matched amphtudes at ki = k2 = = K/3, i.e. in the equilateral limit. 
Defined in this way /nl is essentially the value of the amplitude ^ at a particular point in 
fc-space. Note that we follow the WMAP sign convention here, where positive /nl physically 
corresponds to negative-skewness for the temperature fiuctuations. In the literature expres- 
sions for /nl for a general theory with action (2.1 ) can be found for two cases: (I) For slowly 
varying Cg and to first order in slow- variation parameters [7]. (II) For (potentially) rapidly 
varying Cg when the 2-point correlation function is exactly scale invariant, i.e. Ug = 1 [8J. 
Here we generalize these results for an arbitrary small tilt Ug and to arbitrary order in slow- 
variation parameter e. Hence we will focus on the novel effects on /nl due to the departure 
from scale-invariance and slow-roll. 



An exact, but rather cumbersome, expression for /nl obtained by combining (3.11), (3. 9) 



and (3.10) is given in the appendix A. 2 In the limit of exact scale invariance we recover 



equation 8.4 of |8] for which the following is a fitting formula which illustrates the parametric 
dependence of /nl in this case 

0.164 A + a2 



SI 
NL 



0.27- 



- (0.12 + 0.04/ 



X, 



(3.12) 



where 02 is defined in (3.6). In order to establish a connection between different conventions 



we remind ourselves that fx can be expressed in terms of A/E [55], given that 



/ 



X 



4+^ 



(3.13) 



It follows that |/x| ^ 1 in the presence of cf < 1 entails ^ > 1. From ( |3.12[ ) a good order 



of magnitude estimate for models with a large amount of non-Gaussianity as measured by 
/nl > 1 is 



/nl ~ 0{cf) + 0{ 



(3.14) 



Interestingly, this result corresponds to the case where fluctuations in the scalar field dom- 
inate over those in gravity [59] [72]. Equation (3.14) suggests that deviations from slow-roll 



and breaking of scale invariance alone are not sufficient to generate significant levels of non- 
Gaussianity by themselves. The left graph of figure [T] confirms this, as it shows that, in the 
observationally allowed range, effects from Ug and e alone can yield /nl ~ C^(l) at most. 

However, departures from slow-roll and scale-invariance can have a significant effect for 
non-Gaussianities of detectable size /nl^C>(1)- Note that, considering the parametric de- 
pendence of the amplitude, is fixed in terms of ng,e via (3.6). For a given Ug, an / 



NL 



increasing with larger e therefore corresponds to decreasing whilst /nl increases and vice 
versa. The middle and right graphs in figure [Tj show how considering non-slow-roll e can 
lead to a suppression of /nl- If non-Gaussianity is primarily sourced by Cs ^ 1 (the middle 
graph in[T]), we can see that effects from non-slow-roll can suppress /nl by about a factor 
of 2. Breaking of exact scale invariance similarly has the observationally relevant effect of 
changing the size by up to (i/NL ~ 20. 
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Even more striking is the case when a large amplitude A is primarily due to A/S ^ 
1 depicted in the right graph in figure [T} Then, the exact expression for /nl simplifies 
considerably and we obtain (irrespective of whether Cg is large or small) 

At first sight a good estimate for /nl in this case is /nl ~ 0{fx/c^) ~ C^(^) in agreement 



with (3.14). However, the amplitude is suppressed by a factor of cos^r(3 + 0^2). The 
suppression is so strong that the size of the amphtude can be reduced by more than an order 
of magnitude. In fact, if e is allowed to range all the way up e ~ 0.4, this effect will even 
lead to a change in sign for /nl- Thus slow-roll models with an /nl that violate present 
bounds (see below) can be reconciled with observations when non-slow-roll effects are taken 
into consideration. 

Another interesting limit to consider is that of large speed of sound 3> 1 [Ml EZl EHl EHl 



Q\. Theories with this behavior are addressed in more detail in section 4.3, especially the 
question of how to interpret an apparently "superluminal" speed of sound . Concentrating 
on /nl for the time being, one can expand at first order in — 1 in the limit — )■ 00 to 
obtain 

f^''^'^ ~ 0.28 - 0.04e - (1.19 + 0.08e)(n, - 1), (3.16) 

where we stress that the /^-dependence is suppressed. For a spectral tilt of = 0.96 
this behaves as /nl ~ 0.331 — 0.038e|^ Whilst no /nl ^ 1 is obtained in accordance 



with (3.14), this shows that non-Gaussianities in excess of the levels yielded by slow-roll 
infiation (/^l ~ 0.01) are possible here. 

Observational constraints on /nl depend on which non-Gaussian shape is being probed [751 
[THllTTj, as discussed in the next subsection. Frequently these constraints are quoted for three 
shapes. Firstly, the equilateral shape, peaking in the limit ki = k2 = k2 and going to zero in 
the fiat limit ki + k2 = k^. This roughly corresponds to the shapes following from A^^q^^q^, A, 



^2 



and also approximately (although for ^^^3 there are small differences in the flat limit 
ki + k2 = k^). A second possibility is the local shape, which peaks in the squeezed limit 
ki = k2 ^ ks and goes to zero in the equilateral limit. This is loosely associated with the 
shapes of ^^^2 and ^^(9^)2. Finally, one may consider the orthogonal shape, which peaks in 
the folded configuration 2ki = 2k2 = k^ and goes to zero in the equilateral and squeezed 
limits. This corresponds to the difference between ^^3 and a more exact equilateral shape 
such as A^Qi^Q^ or Ae^. For a list of possible factorizable shape ansatze we refer to [THl [79l |9 



'^Note that whilst a bound /nl > 5 is often quoted |77|[71|, non-Gaussianity from models with lower /nl 
may be detected depending on the shape, if it has a significant amplitude in regions not probed by /nl- As a 
unique shape template is provided by models with Cg 3> 1, deriving constraints for the full particular shape 
may be a promising strategy here and allow detection of low /nl signatures. 

complete basis capable of describing arbitrary shapes in k-space in principle has infinitely many 
"shape-members". However, for the models considered here these three basis shapes will be sufficient. 
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Current WMAP limits on /nl for each of these shapes are P 



fZ = 26 ±240 (95% CL), 
f^l''^ = 32 ±42 (95% CL), 

_ -202 ± 208 (95% CL). (3.17) 



3.3 Shapes 

In the previous section the information contained in A was condensed into a single param- 
eter, /nl- However, there is vastly more information available, encoded in the functional 
dependence of A [751 1151 EZ]- This is uniquely specified by the four parameters e, n^, Cg, /xQ 
In order to disentangle effects from a possible running of A with wavenumber K and the 
shapes themselves, in this section we present the shape of the amplitude at fixed K. In 



interpreting the shapes shown in this section it may be useful to emphasize that from (3.11) 
it follows that a good estimate for the size /nl corresponding to a given plotted shape (we 
normalize one wavevector k and hence plot A{1, ki, / {kik2)) is the value of the plotted 
amplitude in the equilateral limit. 

In the near future the shape and running of A will hopefully become experimentally 
accessible and, if a significant amount of non-Gaussianity is detected, will greatly constrain 
models of structure formation in the early universe. After all, at present these models are 
only observationally constrained to fit or respect upper limits on a number of single value 
parameters (The amplitude of scalar perturbations A, its spectral tilt ns, /nl, the tensor- 
to-scalar ratio r). Fitting a full function would prove far more challenging and provide a 
strong model selection tool. 

The prototypical shape that has been associated with significant non-Gaussianity in 
single-field models is one of the equilateral type [3 El [9l [TOl [TT], where the amplitude peaks 
in the limit ki = k2 = k^. This is in contrast to e.g. multifield models which can also give 
rise to strong local contributions to A (i.e. peaking in the squeezed limit fci ~ A;2 ^ k^) [9] 
via superhorizon interactions local in position space. This is because superhorizon evolution 
is local in space, as different regions are not causally connected with each other. Therefore 
any amount of non-Gaussianity generated after horizon exit is also local in position space 
and hence not local in /c-space. It consequently peaks in the squeezed limit. In single field 
models, on the other hand, ( is frozen in upon horizon exit, so no superhorizon evolution 
takes place. Far inside the horizon modes typically oscillate and their contributions to A 
average out. For known exceptions to this see [m [151 E [IHl HZl HE]- Thus A is almost 
completely sourced by modes exiting the horizon at similar times and hence with similar 
wavelengths. Consequently A peaks in the equilateral limit. 

An interesting question to ask is therefore: Considering deviations from slow-roll and 



^Note that there is an issue of choice here as we have six variables e, es,ns,Cs, fx, -^/S and two constraint 



equations (2.24) and (3.131. For instance specifying e, es,Cs,A/S therefore also fixes the amplitude. 



15 



Figure 2: Here we plot the Non-Gaussian amplitude ^(1, ki, k2)/{kik2) disentangling effects 
from e, Cs and fx- Top row: The effect of varying e. The amplitude is plotted for 
e = 0.01, 0.1, 0.3 respectively from left to right. Cg = 1 and fx = 0.01. Middle row: 
The effect of varying Cs- The amplitude is plotted for Cs = 1, 0.1, 10 respectively from 
left to right, e = 0.01 and fx = 0.01. Bottom row: The effect of varying fx- fx = 
0.01, 100, —100 respectively from left to right, e = 0.01 and = 1. 



exact scale-invariance, how robust is the statement that single field models with sizable non- 
Gaussianity have an equilateral shape? In the slow-roll limit the cubic action (3.1) is domi- 
nated by its first three terms 



low— roll 



M, 2 



I PI 



dtd^x < -a^ 



S I 



+ 2A 



+ 



ae 



6 - 3 + SciXC + ^{e - 26, + 1 - cMdCY + .. 



(3.18) 



This is visible from (3.9) and (3.10) as we can see there that the amplitudes corresponding 
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to the remaining terms in the action are suppressed by powers of e. The leading order 
contributions to A are therefore 



A 

slow— roll 

A^3+A^^2+A^^9Q2. (3.19) 

A breaking of scale-invariance will then introduce a local component into the amplitude. 
This is shown in the left column of Figure |2| where the amplitude is plotted for a slow- 
rolling near-canonical field (i.e. with ~ 1 and fx ~ 0) and spectral index Ug = 0.96. In 
fact it is known that, in the absence of large subhorizon interactions, the very squeezed limit 
is always dominated by a local shape contribution proportional to ris — 1 [7T| 180] 



The middle and bottom row of Figure [2] then show how considering non-canonical fields with 
Cg 1 and/or \fx \ ^ 1 can amplify A. The resulting shapes are predominantly equilateral. 
It is important to note that, whilst the local contribution from breaking of scale-invariance 
can be amplified by non-canonical kinetic terms, the only significant deviation from an 
equilateral shape is in the squeezed limit. The extreme squeezed limit is of course not 
observationally accessible, as it corresponds to considering modes with infinite wavelength. 
If the local contribution is suppressed outside the extreme squeezed limit, as in the case 
plotted in the middle row of Figure [2| it is therefore possible that the observable section 
of the amplitude A is consistent with it being purely equilateral, despite there being local 
contributions. 

What happens once we violate the slow-roll conditions and consider cases with e ~ 
0{1) consistent with constraints derived above? Firstly the remaining interaction terms 
in (3.1),^^g^g^ and A^i, which were previously suppressed by powers of e, now become 



relevant. Even more importantly the remaining terms receive large corrections from e. This 
is most easily shown by considering the case when non-Gaussianities are primarily sourced 
by /x ^ 1- In this case we have 



^(/x»l)~^C3 ~^(^ (e + e.-l)cos^r(3 + a.)^. (3.21) 



As we already saw in the previous section (3.15) this means the amplitude is suppressed by a 
factor of cos ^r(3 + a2). Considering violations of slow-roll can therefore lead to interesting 
new phenomenology, as previously suppressed contributions become important. 

How can one understand the fast-roll suppression of the non-Gaussian amplitude as 



manifest in (3.21) physically? In essence one can separate the effect of considering large 



e, on the amplitude into two categories. On the one hand the interaction Hamiltonian 



that follows directly from (3.1) is modified, since the "coupling constants" for individual 



interaction terms depend on e,es. The interaction terms one can ignore in the slow-roll 



limit (3.18) are all proportional to e, so that their effect linearly increases upon breaking 
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Figure 3: Here we plot the Non-Gaussian amplitude A{1, ki, / {kik2) for fx = 0.01 and 
Cs = 0.1. Plots in the top and bottom rows are for = 1, 0.96 respectively. Plots in the left 
and right columns are for e = 0.001, 0.3 respectively. 

slow-roll. This explains why in the canonical case with ignorable fx (as depicted in the left 
panel of figure [T]) an increase in e leads to an enhanced non- Gaussian signal. For in the 
canonical extreme slow-roll limit one is essentially considering a free field and interaction 
terms "switch on" as larger values of slow-roll parameters are considered. 

The second effect of large e, e^, which lies at the root of the fast-roll suppression of A 



discussed above, is a modification to the propagators (2.21). In other words the functional 
dependence of ( on slow-roll parameters can then become important. The argument is anal- 
ogous to the original one presented for the emergence of an equilateral shape for single field 
models. We recall that this argument invoked modes far inside the horizon oscillating so 
that their contributions to A average out. Considering non-slow- roll propagators can now 



have a similar effect. For the appearance of suppression factors as shown in equations (3.15) 



and (3.21) is essentially due to a "destructive interference" of modes. This means cancella- 
tions between contributions to the amplitude A can occur, since the oscillatory behaviour of 
( is modified. 

Figures |3] and |4] show two concrete case studies to illustrate our results. Figure |3] shows 
how the relative contribution of an orthogonal component is amplified by e and how 7^ 1 




Figure 4: Here we plot the Non-Gaussian amplitude ^(1, ki, k2)/{kik2) for fx = —100 and 
Cs = 0.05. Plots in the top and bottom rows are for = 1, 0.96 respectively. Plots for the 
left and right columns are for e = 0.001, 0.3 respectively. 

induces a local component. The local component away from the extreme squeezed limit is 
also modified by e ~ 0{1). The resulting amplitude gives rise to an intermediate shape with 
both local, equilateral and orthogonal parts. As the bottom right hand graph shows, this 
would be observationally distinguishable from the slow-roll and scale-invariant cases. 

Figure |4] illustrates the e-suppression of otherwise dominant contributions. The top left 
graph shows the amplitude for the given choice of parameters in the case of exact scale- 
invariance and slow-roll. A very large amplitude with /nl ~ 1000 is produced, which violates 
present upper bounds on the level of non-Gaussianity. Deviating away from slow-roll has 
two important effects here. Firstly the fast-roll suppression moves the amplitude back within 
observational constraints. Therefore new regions in the parameter space {cs,/x} open up. 
Secondly a delicate cancellation between terms can bring out contributions from otherwise 
suppressed orthogonal and local shapes even more here than in the example considered 
above. This is shown by the e ~ 0{1) amplitudes in Figure |4] peaking in the folded limit 
2ki 2k2 ~ k^. A shape which cannot be decomposed into local and equilateral shapes 
alone, but requires a strong orthogonal component. 
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Violation of slow-roll and deviation from exact scale invariance therefore naturally lead 
to the generation of intermediate shapes with equilateral, local and orthogonal contributions. 
Hence single field models can produce a richer phenomenology than purely equilateral shape 
type non-Gaussianity. And whilst general statements about limits of the three-point function 
for single field models, such as [TH [80] remain true, such models have a more complex 
fingerprint when considering the full amplitude. 

3.4 The running ting 




Figure 5: Top row: n^Q plotted against e for Ug = 0.96. Green (tzng < 2/3) and yellow 
(riNG < 4/3) regions are those allowed by perturbative constraints on the level of Non- 
Gaussianity assuming |/nl(CMB)| ~ C(IOO) and |/nl(CMB)| ~ C(l) respectively. From left 
to right we plot: {c, = 0.1,/x = 0.01} (left); {c, = IJx = 100} (middle) ; {c, = = 10} 
(right). Bottom left: tt-ng plotted against e in the limit as Cg — t- 0. The black-shaded 
region corresponds to values for n^Q for Ug between 0.9 and 1.1. Almost identical plots are 
found for \fx\ ^ 1. Bottom middle: ting plotted against LogiQ(cs) for Ug = 0.96 and 
e = 0.3,0.1,0.01 from top to bottom respectively, fx = 1- Bottom right: plotted 
against Log]^o(/x) = 0.96 and e = 0.3, 0.1, 0.01 from top to bottom respectively. Cg = 1. 

Having studied the shapes of non-Gaussianity at fixed K in the previous section, we now 
investigate the running of A with scale K [611 |63]. Introducing a dynamical speed of 
sound can in principle lead to a strongly scale-dependent non-Gaussian amplitude A. To see 
why we rewrite Cg, the speed of sound at sound horizon crossing, in terms of the wavenumber 
K (which, we remind ourselves, is defined a.s K = ki + k2 + k^) 

Cg^K^^K (3.22) 
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Different scales K will therefore "see" a different Cg upon crossing the horizon. To quantify 
this difference we follow [61] in defining a spectral index for /nl as 



_ dln|/NL| 

nNG - 1 = , (3.23) 

where we evaluate the running njvc at a fixed point of the amplitude as measured by /nl to 



separate effects from the running and shape. Comparison with (3.9) and (3.10) shows that 
we can always schematically express /nl as /nl = Ci(ns, e, /x)(l — C2{ns,e, fx)c~'^), where 
Ci and C2 are functions of the slow-roll parameters Ci{e,ns, fx) or equivalently Ci(e,es,/ 



X, 



and are given exactly in appendix A. 2 As such we obtain 

In I/nlI = In |Ci| + In |1 - C2cf\. (3.24) 
The non-Gaussian tilt umg is therefore given by 



2 + 2e(-3 + n,) -2n, C2(e,e„/x 
riNG - 1 = 



1 + e cl-C2ie,esJx] 
-2es C2 



(3.25) 



In the small speed of sound limit and for an exactly scale-invariant 2-point function in 
the e > regime the 3-point function acquires a blue tilt [8] 

4e 

nNG - 1 ~ — -. (3.26) 

e + 1 

This is because exact scale-invariance is associated with e<j = — 2e and by requiring e > we 
have prohibited solutions with a phantom equation of state w < —1. However, for slightly 
tilted 2-point functions the behavior of ung is not as simple as we will see. In fact for 
Hs = 0.96 and in the slow-roll limit one obtains a slightly red tilted ung- 

Can we put constraints on the value of the tilt of the 3-point function? In principle levels 
of non-Gaussianity are observable for all scales running from CMB scales {k~^ ~ lO^Mpc) 
down to galactic scales (fc~^ ~ IMpc) [H |62l |8T]. In terms of the scale K this therefore 
corresponds to Kg^^i/ Kqmb — 10'^. For the non-Gaussian tilt this means 

/nl(CMB) ^ 10-3(«NG-i) j^^(Gal) . (3.27) 

Blue/red tilted 3-point functions thus correspond to larger /smaller non-Gaussian amplitudes 
on smaller scales. In order for a perturbative treatment to be applicable the non-Gaussian 
contribution to ( clearly must be much smaller than the Gaussian part /nl ^ — 10^ 
for all observable scales. Present observational constraints on CMB scales give /nl(CMB) < 
(9(100). In case of a scale-invariant or red-tilted 3-point function, where the amplitude is 
largest for large scales, a perturbative treatment is valid all the way down to the smallest 
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scales. However, for blue-tilted 3-point functions one can find interesting constraints on n^G 
and hence on the slow-roll parameters. Specifically we have 

nNG-l<2/3 for |/nl(CMB)| ~ O(IOO) 
riNG-l<4/3 for |/nl(CMB)| ~ (3.28) 

In Figure [5] we show how these bounds translate into constraints on the parameter space 
{?T,s,e, Cs}. The top row shows that the precise bounds on e resulting from tt-tvg are in 
fact model-dependent. The left two graphs show typical plots for large fx and small 
Cs- Very roughly one can see that both scenarios require e < 0.5 to satisfy the weaker 
bound for |/nl(CMB)| ~ 0{1). Whilst this constraint is sufficient for the stronger bound 
|/nl(CMB)| ~ C(IOO) in the case of large /x, this is not the case for a large /nl sourced by 
Cs < 1. There we find e < 0.25 as a requirement. The right graph finally shows that models 
can be tuned to satisfy bounds on n^G for all values in the range < e < 1. Interestingly 
there are also regions with e > 1 where tzng is well-behaved and the constraints are satisfied. 
This may be of importance in non-infiationary model-building. To summarize, we therefore 
find that the bounds on e derived from umg are model-dependent and significantly weakened 
compared to those derived in [8]. 

However, in the limits as fx ^ 00 or very small — )■ these graphs show the same be- 
havior asymptotically. This is shown in the bottom left graph, which furthermore illustrates 
the dependence of bounds on n^. One can read off e < 0.5 and e < 0.25 as approximate 
bounds for |/nl(CMB)| ~ C(l) and |/nl(CMB)| ~ C(IOO) respectively. If we focus on 
the stronger bound for |/nl(CMB)| ~ C(IOO), it is also interesting to note that, whilst the 
constraint for red-tilt with Ug ~ 0.96 is e < 0.25, for a blue tilt we obtain the stronger bound 
e < 0.15 here. 

For simplicity the explicit examples and plots we provide throughout the paper are there- 
fore chosen to satisfy bounds 

e<0.3 , /nl(CMB) < 100. (3.29) 

This ensures constraints from n^G are met for models with slightly red n^. However, we 
stress that choosing larger values of /nl (especially of the equilateral type) at the expense 
of reducing e or fine-tuning models to loosen bounds on e is possible in accordance with the 
expressions derived above. 

The bottom middle graph then reiterates the way in which a large Cg suppresses nNc- 
Strong constraints on the slow-roll parameter e are therefore only obtained for small speeds 
of sound. The graph also nicely shows how we can interpolate between red and blue tilts of 
the three-point correlation function by choosing suitable background variables. Specifically 
the line plotted for e = 0.01 displays a red tilt. The bottom right graph shows the analogous 
plot of riNG vs fx- 

Constraints from n^G can therefore be used to put bounds on e. These bounds can 
become important here, since considering non-slow-roll models naturally gives rise to a large 
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running nj\iG of non-Gaussianity with scale, which are consequently a smoking-gun signature 
of such fast-roll models. Together with measurements of the spectral tilt of the 2-point 
correlation function, measuring hng would allow a direct measurement of e for single field 
models as considered here. Again we reiterate that we use "fast-roll" as a description for an 
evolution which violates the "slow-variation" condition e ^ 1 and not just as a reference to 
properties of the potential ^(0) (see discussion in section [2]). 



4 Concrete model examples 






Figure 6: Here we plot the Non-Gaussian amplitude A{l,ki,k2)/{kik2) for three specific 
models. Ug = 0.96 throughout. Left column: DBI inflation. The amplitude is plotted 
for e = 0.01(top)and e = 0.3 (bottom) respectively. Middle column: Inflation with 
1^1 ^ 1. The amplitude is plotted for e = 0.01(top)and e = 0.3 (bottom) respectively. 
fx = 1000 and Cg = 1. Right column: Bimetric theories with large c^. The amplitude 
is plotted in the limit as — ?■ oo and for e = 0.01(top)and e = 0.3 (bottom) respectively. 
This confirms that A is not sensitive to e for bimetric theories in the large speed of sound 
limit. 



In the previous sections we derived general expressions for the levels of non-Gaussianity 
exhibited by single field models 2T without assuming slow-roll and discussed their associated 
phenomenology. Here we consider three concrete model implementations, writing down 
explicit actions for each model and showing how departure from slow-roll and scale-invariance 
affects their non-Gaussian signatures. 
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For single field models with large non-Gaussianity there are two interesting limits. From 



(3.14) we have 

hL-0{c;') + 0{^), (4.1) 

so that models with ^ 1 and |A/S| ^ 1 will lead to large non-Gaussianities. We also 
considered another case with a smaller, but potentially detectable, non-Gaussian signal /nl ~ 
0{1), when Cg ^ 1. The three concrete examples we provide are therefore examples of 
models with these signatures: (I) DBI inflation [cg ^ 1), (II) Models with A/S ^ 1 and 
(III) Bimetric theories of structure formation [cg ^ 1). 

4.1 Fast-roll DBI inflation 

In DBI inflation [SH] a 3 + 1 dimensional brane moves in warped extra dimensions, giving 



rise to an effective 4D theory with action (2.1), where 



P(X, 0) = -ri(0)v/l-2/(0)X + r 1(0) - y(0) . (4.2) 

/(0) is the so-called warp factor and is constrained to be positive by the signature of the space 
(in this paper we use the convention — h ++). We only consider the 4-D effective field theory 



defined by (4.2). If the full dynamics of a particular higher dimensional implementation are 
considered, far stronger constraints (e.g. from gravitational waves [82]) can be obtained. 
The speed of sound in DBI models is given by 

c, = - 2/(0)X (4.3) 

which is the inverse of the Lorentz boost factor 7, so <^ 1 since the inflaton rolls ultra- 
relativistically. Importantly one also finds 

A = ^(i-a , /r = i-c^, (4.4) 

so that the otherwise typically dominant A^s contribution vanishes, and the amplitude con- 
sequently becomes independent of parameter fx, resulting in: 

- 0{c-% (4.5) 

DBI inflation is the leading example of a single field model with large non-Gaussianity to be 
found in the hterature, and it is associated with an equilateral shape ^ [3, [HI [ID]- Figure |6] 
confirms this for the case of slow roll e <^ 1. As expected the equilateral shape receives a small 
local correction from 7^ 1, which becomes relevant in the squeezed limit. However, just as 



for the more general single field case considered in section 3^, we find that departure from 
slow-roll suppresses the overall amplitude and widens the peaks in the squeezed limit. This 
happens as a consequence of introducing orthogonal and local contributions, thus potentially 
making them observationally accessible. 
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4.2 Fast-roll inflation with ||| > 1 

We now build an explicit action for a model realizing |A/S| ^ 1. First we consider a number 
of constraints that can be placed on the differential properties of P{X,(f)). We can write 
slow-roll parameter e as 

e = —V • (4-6) 

^ XP,x 

In order to satisfy the null energy condition p + p > and 1 < < one must require, 
respectively: 

2XPx>0. , P,xx>0, (4.7) 
as can be seen from the expression for the speed of sound 

p,x Px + 2XPxx 

Observational bounds on parameters such as ns,r,nt therefore only constrain the first two 
derivatives XPx and X'^Pxx [83J. Expressing A and S in terms of derivatives of P with 
respect to the canonical kinetic term X one finds 

^ ^ XP,xx 

SO that in order for |A/S| ^ 1 we require 

\X^P,xxx\ > \P,x\- (4.10) 

Large non-Gaussianities sourced by A/S (or equivalently fx) thus open up the doors for 
constraining P^xxx- 

Adapting the scheme introduced by [83J we can therefore construct a general form for any 
action that gives rise to a large non-Gaussian amplitude A and is consistent with constraints 
on XPx, X'^P^xx and X^P^xxx- We find: 

P(X,0) = g(X,0) + P(Xo,0)-g(Xo,0) 

+ [Px(Xo,0)-g,x(Xo,0)](X-Xo) 

+ ^[Pxx(Xo,0)-g,xx(Xo,0)](X-Xo)' 

+ l[Pxxx{Xo,<P)-q,xxx{Xo,<P)]{X-Xof. (4.11) 
o 

(where q is an arbitrary function of and X) so that higher derivatives along X re- 
main unconstrained. Xq can be fixed by a gauge choice, associated with field redefinitions 
— )■ = g{(j)) (fH3] use the gauge X(Xe) = 1/2 where is the number of e-folds of infla- 
tion). Choosing a value for X at a specific time in this way is also equivalent to choosing a 
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normalization of 0. Having obtained the general action (4.11 ) one can thus write down a the- 
ory with any given value for A/S as discussed in section 3.3, since P^xxx is not constrained 
in any other way than through the value taken by A/S. 

Figure [g] shows that the generic shape generated by ^ 1 is equilateral. In fact much 
more so than in the DBI case, since the n^-dependent local contribution is confined to the 
extreme squeezed limit. As shown in detail in section 3^ deviations from slow-roll strongly 
suppress the amplitude. 



4.3 Bimetric theories 

Bimetric theories [IHl [HI [HI] posit that gravity and matter are minimally coupled to two 
different metrics, related by a single dynamical scalar degree of freedom. Schematically they 
are governed by an action of the form: 

d*x^R[g^,] + / d^x^gC„\g^,,^Matt\. (4.12) 

where g^y and ^^j, are the gravity and matter metrics respectively. The two metrics may 
be seen as independent representations of the local Lorentz group, so that different Lorentz 
transformations must be used to transform measurements made with matter and gravity. 
For this reason no causality paradoxes arise even if the invariant speed of one metric appears 
superluminal with respect to the other metric [l6l [73]. This is in contrast to plain tachyonic 
matter 

Disformal bimetric theories have been considered in the context of cosmological structure 
formation, following from: 

Qt^v = Qt^u - B{(j))df,(j)d^(t) . (4.13) 

It was shown that the minimal such theory (constant, positive B) results in scale-invariant 
curvature perturbations [16], with a distinctive non-Gaussian signature [11]. Departures from 
scale-invariance and their associated non-Gaussian signatures were also examined in [TT] . 
Projecting the scalar field action in the the Einstein (gravity) frame one obtains [^ [TT] an 
action of the DBI type 

S = / Ay=5 (f + Vl + 2W)A-^ - + VW) . (4.14) 

For a solution to the horizon problem, -8(0) = — /(0) > 0, so this is sometimes labelled 
"anti-DBI" (also see [HS]). Obviously a positive B cannot be generated from brane world 
considerations, due to the signature of space-time. But the results derived in the previous 
sections can be easily applied to this class of theories. 

^ A similar argument has been used for single metric, single field inflationary scenarios. If the field serves 
as a time dependent background establishing a preferred coordinate frame, no causality paradoxes come 
about because perturbations only travel superluminally in the preferred frame then [73]. 
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Bimetric models produce a unique non-Gaussian signature. Taking the small tilt {ng—l <^ 
1) and large speed of sound Cs — )■ oo limits, we find [H] 



A 



2K^ 

+ K - 1: 



+0 



1 



(4.15) 



The corresponding expression for /nl is given in |3.16[ Any dependence on e therefore drops 
out and the amplitude is uniquely fixed by the spectral index Us, establishing a consistency 
relationship between and A. This holds great promise for testing bimetric models with 
non-Gaussian constraints, in contrast to the infiationary relation between scalar spectral 
index ris and the tensor-to-scalar ratio r in infiation. Non-Gaussian features in the scalar 
fiuctuations are more likely to be observationally accessible in the near future than tensor 
modes. 

Figure |6] confirms that the amplitude A is independent of e for bimetric models. Note 
that the peaks in the squeezed limit are solely due to breaking of scale-invariance and thus 
are a direct measure of Ug. They may still be observationally elusive, if confined too far into 
the squeezed limit. We emphasize that no significant constraints arise for bimetric theories 
from riNG, as the tilt of the three-point function is strongly suppressed for Cg ^ 1. 



5 Conclusions 

In this paper we have computed the amplitude A, size /nl, shape and running n^c of 
the non-Gaussianity to be found in single field models, without assuming slow-roll or exact 
scale-invariance of the scalar power spectrum. The main results are given in section [3] and 



boil down to equations (3.9) and (3.10), which are general expressions for the non-Gaussian 
amplitude in single field models. All other results, relating to the size, shape and running 
of non-Gaussianities, follow from there. We have shown that observational constraints allow 
significant violations of slow-roll conditions e ^ 1 and es ^ 1 in infiationary single field 
models. As such we derived explicit bounds on slow-roll parameter e in non-scale-invariant 
and non-slow-roll single field scenarios. We found models with e as large as ~ 0.3 satisfying 
all present bounds. 

A much richer non-Gaussian phenomenology can be associated with single field models 
which violate the slow-roll conditions. We established the following results: 

• The size of non-Gaussianity /nl is generically reduced when departing from the slow- 
roll regime. This implies that models with a small speed of sound c^, which violate 
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observational constraints in the slow-roll limit, can still be allowed when considering 
fast-roll scenarios. 



Wherever large non-Gaussianities do arise for single field models, they typically have 
an equilateral shape in the slow- roll limit. We showed that fast-roll models can give 
rise to shapes with larger orthogonal and local contributions (away from the extreme 
squeezed limit). For example, cases where the amplitude A peaks in the folded limit 
2ki = 2^2 = ^3 are constructed explicitly in section |3.3[ Except for large local non- 
Gaussianity near the squeezed limit, which remains a smoking gun of multifield models, 
all the interesting shapes can therefore be generated in single field models. 

Fast-roll models generically give rise to a large running of non-Gaussianity nj^ic with 
scale. However, contrary to models with ng = 1, these can be red-tilted. We predict 



exact values for the running in terms of the model parameters in section 3.4 



• For bimetric models (e > 1) the consistency relation derived in [IT] was confirmed: 
the amplitude ^ is a fixed function of Ug- An equilateral shape was obtained with 
small local corrections in the squeezed limit, proportional to {ug — 1). We also derived 
expressions for /nl, showing explicitly that /nl ~ and positive, and we showed 
that ung ~ for these models. 

• Three explicit model examples of theories that exhibit the phenomenology discussed 
above were constructed in section HI 

In the future we hope to take this work further and investigate a number of issues. Having 
computed general single field expressions, it would be interesting to compute non-slow-roll 
non-Gaussian signatures for further specific models that fall within the single field remit. 
Furthermore, consistency relations between two- and three-point correlation functions exist 
in some limits [HlEOlEIlEHllHIllHnilZainillHS]. In analogy with the bimetric scenario one 
should be able to find constraints uniquely relating observables such as Ug to the full non- 
Gaussian amplitude and shape [93j . This would be a useful extension of results reported in 
the literature regarding limits of the amplitude [TH IHOl EZl [Ml E2]- One could also extend the 
computations presented in this paper to cover slow-roll violations associated with rj ~ 0{1). 
This could help cure problems encountered when embedding cosmological models into e.g. 
supergravity. Finally it will be an interesting task to investigate higher order propagator 



corrections to the three-point function as discussed in section |3.1 

Non-Gaussianity is one of the best tools available for testing theories of primordial struc- 
ture formation. With upcoming experiments promising to measure not just the size, but the 
full functional space mapped out by A, constraints from the shape of non-Gaussianity may 
become even more important in differentiating models. In this paper we have demonstrated 
that single scalar fields can source large non-Gaussianities with a much richer phenomenol- 
ogy than the negative /nl, pure equilateral type non-Gaussianity produced by e.g. DBI 
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inflation. Only time will tell, but while future observations might show that slow-roll single 
field models are ruled out, the simplest explanation could also very well still be a single field 
model. One that violates slow-roll. 
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A Appendix 



A.l Details on the computation of the three-point function 

are 



Here we give some further details showing how the amplitudes presented in section |3.1 
computed explicitly. For further details we refer to [TT] . 

A. 1.1 Mode functions Uk{y) 



In section 2.2 we found that the mode functions for scalar perturbations are given by 



1/2 1/2 
C,s / N Cs 



My) = ^ v,{y) = -^^,Sv^Hi'\-ky)- (A.l) 

aMpiV2e aMpi2'^/^ y e 

We now expand at ^ 1 to obtain an approximate expression for the Hankel function 

Hi'\-ky) = -I ?!£M(zM:![i + ,ky + 0{kyf]e~'^y . (A.2) 

TT 

which gives the following approximate expression for Uk, 

uM - -^ 'IZ^ (^) "^'^ (1 + ^ky)e-'\ (A.3) 

where we have used Y{y f» 3/2) \fT^ 1'^ as scale invariance constrains v — 3/2 ^ 1. One 



can check that modes freeze out in the ?/ — )■ limit by comparing with (2.16) and noting 
that 

^ - Ni/-3/2 



1/2 



{-yy-"' . (A.4) 



Differentiating U}^{y) with respect to y one also finds 



« -> { kh (A.5) 
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A. 1.2 Field redefinitions and the cubic action 



For the final term in the cubic action (3.1), 2/(C)-^|i, we have 



SL 



(A.6) 



/(C) 



— C H CC H 

4c2^ c^H^^ Aam^^ 



idO{dO + d-'{d,d,md,0)] 



1 



[{dOidx) - d-\d.-.d,{d,Cd,x))] , 



(A.7) 



where 5 ^ is the inverse Laplacian. Since f^|i is proportional to the linearized equations of 



motion, it can be absorbed by a field redefinition 

C ^ Cn + /(Cn) . 



(A.^ 



A. 1.3 The individual three-point amplitude computation 

Here we show the computation of ^^^2 in detail. All the other amplitudes are computed 
along the same lines. The corresponding three-point correlator is 



(C(ki)C(k2)C(k3)) 



^^2 = i{2Tl) 5 (ki + k2 + k3)Mfc,(yend)Mfc2(2/end)Mfe3(2/end 

Vend 

X / dy 



^^5^ (e - 3 + 34)ul(y)*^*^ + perm. + c,c. (A.9) 



The subscript "end" means that the quantity has to be evaluated at the end of the struc- 



ture forming (scaling) phase. Making use of (2.16), (A.3) and (A. 5) we take some time- 
independent combinations outside the integral and evaluate them a.t y = ycnd'- 



(C(ki)C(k2)C(k3)) 



s end 



ll/endl 



X 



Vend 



dy{e - 3 + 3c^) — (1 - ikiy)klkly^e'^y + perm. + c.c. . (A. 10) 



— 00+ie 



Note that we have dropped a factor of 11^(1 + ikjyendjG'^^^"""^ as this will be negligibly 
small in the limit \ky\ <^ 1 where the truncated Hankel function expansion (A.2) is valid. 
Evaluating the integral and re-expressing variables in terms of quantities calculated at sound 
horizon crossing (when, by convention, y = K~^) we obtain 



(C(ki)C(k2)C(k3)) 



{27rf6%ki + k2 



HHe + e. 



ir2 



(n,A;3^3-2. 



~1{ 



(e — 3) cos — — r(l + a2) 



l + (l + a2 



h 

K 



16Mp;\2 c4 

+ 3c^ cos — — r(i 



Oil, 



^ K 



+sym. 
(A.ll) 
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where, as defined before in (3.6), 



ai = 3 - 2z/ = ■ ; 02 = , -■ 



(A.12) 



A. 1.4 A useful integral 



Here we report a results from [8] for evaluating integrals of the type (A. 10). These integrals 
are typically of the form 



C 



Vend 



dy 



Vend 



-lyYe 



n iKy 



(A.13) 



For 7 + 72 > —2 the imaginary part of (A.13) is convergent as ?/end — > 0. Note especially the 
ie term which is necessary in order to regularize the integral. We can now approximately 
extend the upper limit of integration to 0, which amounts to neglecting terms of higher 
order in (/c|?/end|) consistent with the Hankel function expansion (A.2) we adopted earlier. 
One finds 

(A. 14) 



77r 



ImC = -(i^li/endl)-^ COS ^r(l + 7 + n)ir-"-i 



A.2 General expressions for /nl and n^vc 



Here we give a general expression for /nl as derived from the full amplitudes (3.9 ) and (3.10 ). 
We remember that /nl was defined as 



/l 



NL 



30- 



(A. 15) 



Substituting in the correct amplitude one obtains 



r general 

/nl ~ 

g . 2-3-ns2-2— Sna 

c2(l + e)(n,-2) 



36c^(l + e)(-27 + n,(13 + 7n,))rK]sin 



UsTl 



+ 



16c^(l + e)2r[2 + n,]sin 



L 2 



8e (21 - 22n, + Sn,^) - 4(33 + n,(5n, - 31))) T 

4 — ie — Us + Sen 



(-9 {9e%ns - 2){n, - 1) - 36^(38 + 5{n, - 7)n,) + 
(3e-l)(n,-2)" 



16(l + e)^(/x-l)r 



1 + e 



sm 



1 + e 
3ns) + ns)Tr 



+ 



2(1 + e) 



We also found that /nl can be expressed as 



ln|/ 



NL 



In ICil + In 11 — C2C 



so that the non-Gaussian tilt n^^c "was given by 



TlATG - 1 



-2e. 



e,Ar e-\cl-C2' 
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(A.16) 



(A.17) 



(A. 18) 



where Ci — Ci{ns, e, fx)- These functions are 



Ci{ns,e,fx) 



5(-243 + n,(121 + 67n, + 46(1 + n,)))r[n,]sin [^] 

2l+"^32+3n.(^^ - 2) 



(A.19) 



C2{ns,e,fx) 



{l + e){n -2) ^^''^""^ ~ ^^^""^ ~ ~ ^''^^^ + ^^""^ ~ ^^""^^^ 
8e (21 - 22n, + Sn^^) - 4(33 + n,(5n, - 31))) T 




ifiM ^ ^2/. r4-4e-ns + 3enJ\ . 
16(1 + e) (/x - l)r ^-p- 1 sm 



(A.20) 
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